Chiral Symmetry Breaking and the Dirac Spectrum at Nonzero Chemical Potential 
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The relation between the spectral density of the QCD Dirac operator at nonzero baryon chemical 
potential and the chiral condensate is investigated. We use the analytical result for the eigenvalue 
density in the microscopic regime which shows oscillations with a period that scales as 1/V and 
an amplitude that diverges exponentially with the volume V = L"^. We find that the discontinuity 
of the chiral condensate is due to the whole oscillating region rather than to an accumulation of 
eigenvalues at the origin. These results also extend beyond the microscopic regime to chemical 
potentials p ~ 1/L. 



Introduction. One of the salient features of QCD at low 
energy is the spontaneous breaking of chiral symmetry 
characterized by a discontinuity of the chiral condensate. 
More than two decades ago it was realized by Banks and 
Casher Q that the discontinuity of the chiral conden- 
sate at zero quark mass is proportional to the eigenvalue 
density of the QCD Dirac operator. This relation estab- 
lishes that the eigenvalue spectrum of the anti-Hermitian 
Dirac operator becomes dense at the origin of the imagi- 
nary axis in the thermodynamic limit. The Banks-Casher 
relation is of substantial practical value for nonperturba- 
tive numerical studies of QCD. It allows one to extract 
the chiral condensate directly from the spectral density. 

At nonzero baryon chemical potential, ^, the Euclidean 
Dirac operator, D = DriJ-q + M70j is non-Hermitian so 
that the support of its spectrum is a 2 dimensional do- 
main in the complex plane. In this case we will show that 
the discontinuity of the chiral condensate at zero mass is 
not due to the accumulation of eigenvalues near zero. 

The QCD partition function at zero temperature does 
not depend on the baryon chemical potential. 



ZNf{m;fi) = ZNf{m;^jL = 0), 



for 



< Aic, (1) 



where is the mass per unit quark number of the lightest 
excitation with nonzero quark number. Here and below 
we only consider the case with Nj quark flavors with 
equal mass m. Therefore, the chiral condensate given by 



1 



NfV 



(2) 



remains equal to its value at /i = for /i < /i^. Our aim 
is to understand this behavior from the spectrum of the 
Dirac operator. 

We will consider two types of gauge field averages, 
quenched averages where the determinant of the Dirac 
operator is not included in the average and unquenched 
averages which include the fermion determinant. The 
quenched spectral density can be expressed as 



PQ{x,y;f^) = Cy^S'^ix + iy - zk)), 



k 



(3) 



where the eigenvalues of the Dirac operator are given by 
Zfc, and the brackets denote the average over the Yang- 
Mills action. The eigenvalue density of full QCD includes 
the fermion determinant in the average 



PNfix,y,m;fi) 



(Efe ^Hx + W- zk) det^^ {D + m)) 



(det^^p + m)) 



(4) 



Strictly speaking, since p^j is in general complex, this 
is not a density. Due to chiral symmetry the eigen- 
values occur in pairs ±Zfe so that pNf{x,y,m;n) = 
PNf {-~x, —y, m; p). A second reflection symmetry which 
holds only after averaging over the gauge field configura- 
tions is that (a;, y, m; p) = p^j (x, —y, m; p). Because 
the fermion determinant vanishes for — ±771 we ex- 
pect that pNf{x = irn,y = 0,m;p) = 0. The chiral 
condensate in the chiral limit can be expressed as 



S = lim lim — 

m^O V^oc V 



dx dy 



PNfix,y,m; p) 
X + iy + m 



(5) 



At zero chemical potential the above quantity is known 
to be proportional to the density of the imaginary 
eigenvalues at zero 1|. This happens because then 
Pn f{x,y,m; p = 0) oc 5{x). When p ^ Q the eigenval- 
ues spread into the complex plane and this argument no 
longer holds. In this case we will show that there is an ex- 
tended region of the eigenvalue density that contributes 
to the chiral condensate. 

For simplicity, here we will show how the condensate 
arises from the microscopic limit of the spectral density 
which is believed to be universal. We have also verified 
that a similar mechanism occurs even for larger p 
provided that the contributing eigenvalues are still in the 
universal region. The microscopic limit of the spectral 
density is defined as 



PNi{x,y,rn\p) 

1 . X y rh 

' P^f^s^' vt7' vt7' 



(6) 
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^sv sy' F^Vv 



In this limit, x — xT,V , y — y'SV, rh — mY,V and p = 
pFttVV are kept fixed with E given by ISJ and iv the 
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pion decay constant. The expression for the condensate 
(Isl) in this hmit becomes 



hm E 



dx dy 



PNf{S:,y,m;fi) 



ly 



(7) 



The microscopic hmit of the spectral density at nonzero 
chemical potential was recently calculated both for the 
quenched case and the unquenched case For 
a nonzero number of flavors it was found that the 
eigenvalue density for mS < Ifj^F^ is a strongly oscil- 
lating complex function. The oscillations cover a region 
of the complex eigenvalue plane and, as we will see be- 
low, the entire region contributes to the integral in Q. 
This constitutes a new mechanism where a discontinu- 
ity of the chiral condensate in the complex mass plane 
is obtained from an oscillating eigenvalue density in the 
complex plane. This mechanism does not rely on the spe- 
cific form of the eigenvalue density as is demonstrated in 
the simple example below. 

The lack of Hermiticity properties of the Dirac op- 
erator at nonzero chemical potential is a direct conse- 
quence of the imbalance between quarks and anti-quarks 
imposed in order to induce a nonzero baryon density. 
Because of the phase of the fermion determinant, proba- 
bilistic methods are no longer effective in the analysis of 
the partition function. This is known as the sign prob- 
lem. Although progress has been made in some areas we 
believe that because of its physical origin, a paradigm 
shift will be necessary to develop viable probabilistic al- 
gorithms for this problem. Because of this, it is our opin- 
ion that it is particularly important to improve our an- 
alytical understanding of chiral symmetry breaking for 
QCD at nonzero baryon density. 

Euclidean QCD at finite baryon density is not the only 
system without Hermiticity properties that has received 
much attention recently. We mention the distribution of 
the poles of S— matrices which are given by the eigen- 
values of a non-Hermitian operator 0, 0| , the Hatano- 
Nelson model Q (a random potential together with a 
nonzero imaginary vector potential) , and the description 
of Laplacian growth in terms of the spectrum of non- 
Hermitian random matrices fioj . The essential difference 
from QCD is that in these problems the determinant of 
the operator only enters in the generating function of the 
resolvent. We will see that the additional determinant in 
QCD completely changes the character of the theory. 

Example. As an example to illustrate our point, let us 
consider the eigenvalue density (sg is the sign function) 



PE^{x,y,m;fi) 



(8) 



X 1 _ gisg(«)(l*l+2A")S/4A' + (|41-l'5i|)(l*l+2A')/4A^ 

This eigenvalue density has the same reflection symme- 
tries as the eigenvalue density of the QCD Dirac operator 




. 00, >/ 

-100^ — 



Re[/5jv^=i(a;, y,m;/i)] x 

FIG. 1: Top view of the real part of the eigenvalue density 
for one flavor with m = 60 and /i = 8 in half of the complex 
eigenvalue plane. The oscillations are cut off to illustrate the 
eigenvalue repulsion at x = y = and the drop off for x > 2fi^. 

and has the property that it vanishes at the point where 
the fermion determinant is zero. The integral in |2J can 
be evaluated analytically by means of a complex contour 
integral in y resulting in 



I]Ex = sg(TO)S-h-e-l"l(e- 
m 



-\m\ 



I) 



(9) 



which, for large m, approaches sg(m)S. What we have 
learned from this example is that a discontinuity in the 
chiral condensate can be obtained from an oscillating 
spectral density rather than from eigenvalues localized 
on the imaginary axis. We will show next that the same 
mechanism is at work for QCD at ^ =/= 0. 

The microscopic spectral density. The input for our 
calculation of the chiral condensate is the microscopic 
spectral density derived for any number of flavors in 
5] . Without loss of generality we will consider the case 
Nf = 1 and topological charge equal to zero for which the 
equations are less extensive. The microscopic eigenvalue 
density can be decomposed as 0, |^ 

PNf=i{x,y,rh;fi) = pQ{x,y;fi) ~ pu {x , y , m; fi) , (10) 

with {z — X + iy) 

Pa(£,y>;/i) = ^e-(^^+^*^)/(sA^) (11) 



x^o(^) 



dtte-'^'^^*^ Io{z*t)Io{mt). 



The first term in (I10|l is the quenched eigenvalue density 
given by pq{x, y; p) = pu{.x, y, x+iy; fi). As expected, 
the microscopic spectral density vanishes at z ±to. A 
plot of the real part of the eigenvalue density for to = 60 
and p = Sis shown in figure^ Notice that the oscillatory 
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region extends from the mass pole at z = rh and toward 
the boundary of the support of the spectrum. 

The oscillations appear as the microscopic variables 
become large, i.e. as the thermodynamic limit is ap- 
proached. In this region an asymptotic formula for the 
eigenvalue density is accurate and will be used in order 
to analyze the role of the oscillations for chiral symmetry 
breaking. We first derive the asymptotic formula for the 
eigenvalue density and then evaluate the chiral conden- 
sate from (0). For /t^ > 1 and (x -I- rTi)/(4/i^) < 1, the 
integral in Hll|l over t is very well approximated by 



dtte~^^'^*^Iaiz*t)Ia{mt) 



(12) 



Furthermore, we are interested in the approach to 
the thermodynamic limit where \rh\ ^ 1, \z\ 3> 1, 
|fp/(4/i2) > 1 and \rhz\/{Aii^) > 1. This justifies the re- 
placement of the Bessel functions by their leading order 
asymptotic expansion including the Stokes terms. We 
obtain the following asymptotic result for the difference 
between the quenched and the unquenched eigenvalue 
density 



(13) 



J_p-[y+«Bg(£)(Ix| + Im|-4/i=)]V(8/i")-(|£t-2A")V(2A') 



This expression has the reflection symmetries discussed 
below The asymptotic expansion of the quenched 
part of the spectral density is simply given by 



/5g(i,y; A) 



(14) 



For the argument presented below it is important that 
also the asymptotic expansion of the spectral density van- 
ishes at X + iy = ±rh. 

The chiral condensate. As explained in the introduction, 
the chiral condensate does not depend on the baryon 
chemical potential. Hence, in the microscopic limit it 
is known that 0, 0| (momentarily we use the original 
variables to emphasize the volume dependence) 



(15) 



which is discontinuous, TiNf=i{rn) = sg(m)E, in the ther- 
modynamic limit at fixed quark mass. The question we 
wish to answer is how the oscillatory spectral density 
conspires into a fj, independent condensate. We stress 
that this is not just a challenging mathematical problem; 
understanding which parts of the eigenvalue density con- 
tributes to the chiral condensate will give direct insight 
in the physical consequences of the sign problem. 



We now derive the chiral condensate from lO using the 
asymptotic form of the microscopic eigenvalue density. 
We first consider the integral over y. The contribution 
from the quenched part of the spectral density 114|) is 
given by 



Airfj,^ J_ 



dy 



1 



X + iy + m 



sg{x + m)j^. 



(16) 



The contribution from pjj in Q is evaluated by a saddle 
point approximation. The contour in the complex y plane 
is deformed into a contour from — oo to oo over the saddle 
point at y = isg(a;)(4/i^ — — |m|) and, if the contour 
has crossed the pole, an integral around the pole at y = 
i{x + rh). The saddle point contribution is exponentially 
suppressed for |a;| < 2fi^ leaving only the integral around 
the pole. For m > (m < 0) the pole contribution for 
X > (i < 0) is exponentially suppressed. We obtain 



1 



47r/t2 J 



dye 



-[a+isg(i;)(lx| + lm|-4A")]V8A'-(|Sl-2A')V2A' 



1 



X + iy + rh. 



(17) 



1 



-^^[9{rh)0{—x — rh) — 9{~rh)9{x + rh)], 



where we have used that the exponent vanishes at the 
pole. For x > 2//^ the eigenvalue density is zero so it is 
now trivial to do the integration over x to get 



|,|_^_^dx[isg(£-K,n) (18) 

+9{m)6{-x -rh) - 9{-m)0{x + rh)] 
sg(m)S. 



This result agrees with IjlSII for | rfi | ^ 1 where the asymp- 
totic expansion of I]jvy,=i(m) is valid. Using the exact 
microscopic spectral density we would have recovered the 
mass dependence of (|15|l . 

We also emphasize that a finite result for the chiral 
condensate is not obtained due to a cancellation of the 
pole and a zero of the fermion determinant. The pole 
term gives a finite contribution for each of the two terms 
in (|10|l which do not vanish at z = rh. We have checked 
numerically that the same mechanism results in a discon- 
tinuity of the chiral condensate for more than one fiavor. 

The contribution from the unquenched part of the 
eigenvalue density to the chiral condensate is dominated 
by the pole term because the exponential in (|17ll sup- 
presses the integrand at the saddle point in the complex 
y-plane. The simple result (|17(l which implies a nonzero 
chiral condensate in the chiral limit is thus directly re- 
lated to the complex phase of the spectral density. The 
oscillating exponential has to suppress terms that diverge 
exponentially with the volume which is achieved by oscil- 
lations in the y-direction with a period that scales as the 
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inverse of the volume. In contrast, for quenched QCD, 
where the eigenvalue density is real and positive, the chi- 
ral condensate vanishes in the chiral limit provided that 
/i 7^ as follows immediately from H16|) . Instead a di- 
quark condensate forms and breaks chiral symmetry just 
like in phase quenched theories (see e.g. Chiral 
symmetry breaking in phase quenched and full QCD oc- 
cur by means of two different mechanisms, and the oscil- 
lations of the eigenvalue density due to the sign problem 
in the unqucnchcd case distinguish between the two. 

The resolvent. One can easily convince oneself that for 
z ^ rh the exponent in the integrand of the resolvent 

/V ^* ^ V [ -^2 - PNf (Re(u) , Im(M) , to; /i) 
Gnj [z, z ,m; 11) = h / a u — '' ^ " ' ' 

results in expressions that diverge in the thermodynamic 
limit. Indeed, this was concluded from earlier Random 
Matrix calculations of the resolvent Using the mi- 
croscopic spectral density one finds good numerical 
agreement with the results obtained in |l6| . 

Phase transitions in generating functions. Our results 
can be understood in terms of phase transitions for gener- 
ating functions for the spectral density. Using the replica 
trick Idllil^ 

PNf{x,y,m;fi) = lim — 9^.9^ log Z^v,. ,„(m, z, z*; /x) 

we are naturally lead to the generating functionals 

ZNf,nim,z,z*;fi) = (det(i:>-|-m)^^|det(i:) + z)p") 

for the eigenvalue density. The presence of conjugate 
quarks in the generating function induces a coupling to 
the chemical potential in the low-energy eff'ective theory 
which is completely fixed by the pattern of chiral sym- 
metry breaking and determines the eigenvalue density 
uniquely in the microscopic regime. In order to calcu- 
late the microscopic spectral density this way, one has to 
employ powerful integrability relations that exist for the 
effective partition function pLslligj . 

Let us compare the phases of the generating function 
for n = 1 to the regions in the complex plane charac- 
terized by a different behavior of the eigenvalue density 
for Nf — 1. For real z the phase structure follows from 
|2fll |. The extension for complex z can be obtained from 
the asymptotics of the microscopic generating functions 
given in [6j. For |Re(z)| > 2fP the generating functional 
is in the normal state while for |Re(i)| < 2fi^ two Bose 
condensates are separated by a first order phase tran- 
sition. This phase transition occurs exactly where the 
oscillating region of the eigenvalue density begins, while 
the transition to the normal phase corresponds to the 
drop off of the eigenvalue density for |Re(z)| > 2jj,^. 

Conclusions. For QCD at nonzero chemical potential the 
chiral condensate is not dominated by the contribution 



from the smallest eigenvalues. On the contrary, we have 
found that the contributions from strips parallel to the 
imaginary eigenvalue axis do not depend on the real part 
of the eigenvalue as long the eigenvalue is inside the sup- 
port of the spectrum. This result arises from integrating 
a spectral density that oscillates with a period of 1/{T,V) 
and an amplitude that diverges exponentially with the 
volume. Although we have shown only results using the 
microscopic eigenvalue density, we have checked |3| that 
our arguments apply up to /i ~ 1/L (for V = L'^). In con- 
clusion, we have uncovered a novel mechanism of chiral 
symmetry breaking at nonzero chemical potential where 
an oscillatory spectral density results in a discontinuity 
of the chiral condensate in the complex mass plane. 
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